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We analyze the long-time quantum dynamics of degener- 
ate parametric down-conversion from an initial sub-harmonic 
vacuum (spontaenous down-conversion). Standard lineariza- 
tion of the Heisenberg equations of motions fails in this case, 
since it is based on an expansion around an unstable clas- 
sical solution and neglects pump depletion. Introducing a 
mean-field approximation we find a periodic exchange of en- 
ergy between the pump and subharmonic mode goverened 
by an anharmonic pendulum equation. From this equation 
the optimum interaction time or crystal length for maximum 
conversion can be determined. A numerical integration of 
the 2-mode Schrodinger equation using a dynamically opti- 
mized basis of displaced and squeezed number states verifies 
the characteristic times predicted by the mean-field approx- 
imation. In contrast to semiclassical and mean-field predic- 
tions it is found that quantum fluctuations of the pump mode 
lead to a substantial limitation of the efficiency of parametric 
down-conversion. 



I. INTRODUCTION 

Owing to its relative simplicity but yet richness, the 
process of parametric down-conversion is one of the most 
intensively studied in quantum optics Here pho- 

tons of a coherent pump field are transformed into pairs 
of signal and idler photons (^J^ which can display non- 
classical quantum correlations || or perfect squeezing in 
the case of degeneracy. We here restrict ourselves to the 
latter situation, where both down-converted photons are 
emitted into the same radiation mode. A standard ap- 
proach to analyze the quantum fluctuation in nonlinear 
optical system is to assume small fluctuations around 
the classical solutions, i.e. to linearize the Heisenberg 
equations of motion . The linearization approximation 
fails however in the case of a vacuum input of the sub- 
harmonic mode, since it neglects pump depletion and is 
thus only valid for an infinite input intensity of the pump 
field. Thus linearization can neither be used to study the 
effect of finite system size, i.e. finite pump intensity nor 
the long-time dynamics of the parametric process. 

Using a short-time perturbation expansion, Crouch 
and Braunstein analyzed the leading order corrections 
to the maximum degree of squeezing due to finite pump 
intensities 0. Here we are interested in the long-time 
behaviour of parametric down-conversion. In particu- 



lar we aim to determine the optimum interaction time 
(propagation length in the crystal) for maximum down- 
conversion and the maximum efficiency of this process. 
In the case of a vacuum input of the sub-harmonic mode, 
both quantities are goverened by quantum effects. We 
find that in contrast to the classical predictions, these 
quantum effects limit the maximum conversion efficiency 
from a pump photon into two sub-harmonic photons to 
a value much less than unity. This limitation could be of 
importance for applications in quantum communication 
and cryptography on the single photon level. 



II. MODEL, CLASSICAL DYNAMICS AND 
LINEARISATION 

In order to describe stationary parametric conversion 
of travelling-wave pump radiation into travelling-wave 
sub-harmonic radiation we introduce a moving coordi- 
nate system. Ignoring transversal degrees of freedom we 
find the following Heisenberg equations of motion 



— <zi = Ka 2 a\, 

d K* 2 

dt a2 = ~— 



(1) 
(2) 



ai and a 2 are the bosonic mode operators of the sub- 
hamonic and pump fields respectively. K describes the 
strength of the nonlinear process. It is proportional to 
the nonlinear susceptibility x °f the crystal and the 
inverse of the beam diameter. The time evolution in the 
moving frame corresponds to a spatial evolution in the 
lab frame and the fields at t = are the input fields. Due 
to the phase symmetry of the equations 
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we may choose K and the initial amplitude of the pump 
field (a2(t = 0)) real. The equations of motion (QJ|) 
obey the Manley-Rowe relation HQ, which states that 
the total energy of the free (!) system is conserved. 



— (ala 1 }+2—(ala 2 ) = 0. 



(3) 



Even though Eqs.(|l]) and (|^) seem simple, the non- 
linearity prevents an analytic solution of the quantum 
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problem. Therefore approximations are necessary. A fre- 
quently used approximation is the linearisation around 
the classical solutions. In order to discuss the validity 
of this approximation, let us first consider the classical 
problem, where the Bose operators a\ and a 2 are replaced 
by c- numbers a\ and a-i- 



—at = Ka 2 a\, 
d K 2 

dt a2 = ~Y a i- 



(4) 
(5) 



One clearly sees, that for vanishing sub-harmonic input, 
i.e. ai(0) = 0, both amplitudes remain constant. This 
solution is linearly unstable and any fluctuation will be 
exponentially amplified. The time evolution critically de- 
pends on the amplitude and phase of an initial classical 
fluctuation. Thus a classical calculation cannot deter- 
mine the optimum interaction time (or crystal length) 
for maximum conversion. 

In the standard linearization approach, the pump- 
mode operator is replaced by its classical input ampli- 
tude. This turns the quantum problem into a linear 
one, which can immediately be solved. One finds that 
the time-evolution operator of the sub-harmonic mode is 
given by 

iMt) = s[ v (t)], (6) 

where S is the so-called squeezing operator |TlJ 

SM=exp{|al 2 -^a?} (7) 
with a squeezing parameter that grows linear with time 



r)(t) = Ka 2 t. 



(8) 



Since K and a 2 have been choosen real, the time evo- 
lution will lead to a squeezing of the fluctuations of the 
out-of-phase component of the sub- harmonic mode pi, 
(ai = x\ + ipi) below the standard vacuum limit. The 
quantum noise of pi monotonously decreases with time 
and simultaneously the quantum noise of x\ increases. 
The increase of the fluctuations in the in-phase compo- 
nent x\ is associated with a steady increase of the sub- 
harmonic photon number 



(a\ai) — sinh 2 rj(t) — sinh 2 Ka 2 t. 



(9) 



This result violates the Manley Rowe relations (g) and 
indicates the breakdown of the linearization for larger 
times. The growing fluctuations of the (anti-squeezed 
component of the) sub-harmonic mode can at some point 
not assumed to be small anymore. They will lead to a 
decrease (depletion) of the pump-mode amplitude and to 
fluctuations in this mode. 



III. MEAN-FIELD APPROXIMATION AND 
OPTIMUM INTERACTION TIME 

As noted above a linearization of the Heisenberg equa- 
tions of motion cannot be used to study the long-time 
behaviour of spontaneous (vacuum input) parametric 
down-conversion. The quantum fluctuations of the sub- 
harmonic mode and their backaction onto the pump 
mode are essential and need to be taken into account. 
We may however replace the pump-mode amplitude by 
its average value, which amounts to a mean-field approxi- 
mation p2| . With this we obtain the equations of motion 



— ai = K (a 2 )a\, 
d / \ K * 1 2\ 



(10) 
(11) 



Thus we have transformed the original set of nonlinear 
operator equations into a linear operator equation plus a 
nonlinear classical one. One easily verifies that equations 
(|l0|) and (fTO) obey the Manley-Rowe relation. 



^(ajoi) = K{al)(a\) + c.c. = -2^<4<z 2 >. (12) 

The mean-field equations correpond to a time-evolution 
operator 



U mf (i) = D 2 [/?(*)] Si [v(t)], 



(13) 



that consists of a coherent displacement operator for 
the pump mode and a squeezing operator for the sub- 
harmonic mode. 

D(a) = expjcW — a*a}, 

S( 7? )=exp{|at 2 -^a 2 }. 

Thus the interaction leads to a shift of the coherent am- 
plitude of the pump mode by the amount 



0{t) = ~K j&t'(a\{t')). 



(14) 



At the same time the sub-harmonic mode is squeezed by 

rt 



V (t)=K dt' (a 2 (t')). 
Jo 



(15) 



In contrast to the linearisation, the squeezing parame- 
ter does not increase indefinitely, since the pump mode 
amplitude decreases, characterized by the displacement 
parameter (3. (3{t) and rj(t) are not independent. From 
the mean-field equations we find 



(16) 



If we know rj(t) we can immediately obtain the ampli- 
tude of the (classical) pump mode from Eq. ( |l5| ) . On the 
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other hand we find the following coupled equations for 
the sub-harmonic photon number and correlation func- 
tion 



dt 



(ajai) = 2t) (aiai), 
(axat) = j) (2(a\ai) + l\ 



which have the solutions 



(a\ai) = sinh 2 77, 
(aiai) — — sinh 2ij. 



(17) 
(18) 



Thus the knowledge of 77 is sufficient to determine all 
relevant quantities. From ( |l4]) and (|l6|) we find 77 = 
— (K 2 /2)(a1) and thus the dynamics of the squeezing pa- 
rameters is goverened by an anharmonic pendulum equa- 
tion. 



i){t) = --Jf 2 sinh2r;(7;), 



(19) 



with the initial conditions 



17(0) = 0, 
f l (0)=K{a 2 (0)}. 

The anharmonic pendulum equation with the given ini- 
tial conditions is equivalent to the integrated Manley- 
Rowc relation 



• 2 • 1 2 
—ri] + sinh 77 
K z 



2n 2 ! = 2|(a 2 (0))| 2 



(20) 



This suggests a mechnical analogue. If 77 is interpreted as 
the spatial coordinate of a classical particle moving in one 
dimension, the first term in Eq.(20) represents its kinetic 
and the second its potential energy. In the chosen units 
the kinetic energy is then twice the pump-mode photon 
number and the potential energy the photon number of 
the sub-harmonic mode. 




K{a 2 (Q))t 

FIG. 1. Squeezing parameter as function of scaled time for 
coherent pump with input intensity (02(0)) 2 = n^. 



Fig. 1 shows the squeezing parameter as function of the 
scaled time K (ci2(0))i for different initial photon num- 
bers 7I2 := I (02(0)) | 2 . The squeezing parameter reaches 
a maximum value and there is an optimum interaction 
time or crystal length for maximum squeezing. The cor- 
responding optimum time is a quarter of the oscillation 
period in the anharmonic potential. 



n 2 ~~ 2 sinh 2 y 



(21) 



where sinh y n 



2ri2- This results agrees with that 



of the short-time perturbation expansion by Crouch 
and Braunstein ||. A comparision with the Crouch- 
Braunstein result shows however that the maximum 
amount of noise reduction found in mean-field approx- 
imation 
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32n 2 ) 



(22) 



is too small. The mean-field approach neglects the fluc- 
tuations of the pump mode, in particular its phase noise. 
When this is taken into account the minimum fluctua- 
tions are only (Apf) min = 1/8^1 §@- 

Maximum conversion of pump into sub-harmonic pho- 
tons is achieved when 77 = 0, i.e. at the turning points of 
the classical pendulum motion. Thus the optimum con- 
version time T conv or equivalently the optimum crystal 
length is determined by 



KT — 



dy 



1 



In 



i sinh 2 y 



(23) 



which is twice the time of maximum squeezing. Fig. 2 
shows the scaled photon numbers of the pump and sub- 
harmonic mode as a function of time. 
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K{a 2 {0))t 

FIG. 2. Scaled mean photon number of pump (dashed) 
and sub-harmonic mode (line) as function of scaled time for 
(o 2 (0)} = V20 
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Since the mean-field approach takes into account the 
energy transfer from the pump mode into sub-harmonic 
fluctuations, it correctly describes the oscillatory energy 
exchange in parametric down-conversion from an initial 
sub-harmonic vacuum. This is in contrast to the classi- 
cal or linearization approximation. The mean-field ap- 
proximation also allows to determine the optimum inter- 
action time for large squeezing or best down-conversion, 
Eqs.(^l|,^3|). The maximum conversion efficiency is unity. 

The underlying assumption of the mean-field approach 
is a quasi-classical description of the pump field. This as- 
sumption becomes however questionable at the point of 
total energy conversion and thus the maximum conver- 
sion efficiency obtained in mean-field approximation may 
not be correct. To calculate this quantity and to discuss 
the influence of quantum fluctuation in particular of the 
pump mode we shall numerically integrate the two-mode 
Schrodinger equation in the next section. 

IV. NUMERICAL INTEGRATION OF 
TWO-MODE SCHRODINGER EQUATION AND 
QUANTUM LIMIT TO THE CONVERSION 
EFFICIENCY 

A direct numerical integration of the Schrodinger or Li- 
ouvillc equation is not a straight forward task for multi- 
mode problems. Unless the interacting modes contain 
only very few photons, the standard Fock-basis expansion 
requires the use of a large basis set. For the present prob- 
lem a large basis set is required in both modes since dur- 
ing the interaction all or almost all photons of the pump 
mode are converted into sub-harmonic photons and vice 
versa. 

To avoid the large-memory requirement of a simple 
Fock space expansion one may think of choosing a mod- 
ified basis adapted to the problem. For example in the 
initial phase of the process the pump mode is in a coher- 
ent state |a°). Its photon number distribution is Poisso- 
nian and thus the required number of basis states is of 
the order of la^l, which can be large. On the other hand 
one can displace the number state basis with the unitary 
transformation D(a) introducing the states 

\a,n) = D{a)\n) (24) 

which form a complete set. Clearly at t — only a single 
state is needed to describe the pump mode if a = a®. 
As known from the mean-field approach, the coherent 
amplitude of the pump mode decreases during the inter- 
action and the basis set ( |24| ) would soon become inef- 
fective. Thus the parameter a needs to be dynamically 
adapted, a — > a(t). This is easy to implement in a nu- 
merical algorithm that solves the differential equation. In 
each time step the expansion coefficients are calculated 
in an adapted basis which uses parameters obtained in 



the previous time step. These coefficients are then used 
to update the basis and so on. 

If there is no initial symmetry-breaking the coherent 
amplitude of the subharmonic mode remains zero at all 
times. Thus a dynamically adapted coherent displace- 
ment of the sub-harmonic basis states is not useful. How- 
ever we have seen in the previous section that the time 
evolution of this mode is approximately described by a 
dynamical squeezing S(?y), see Eq.([l3"|). Therefore we ex- 
pand the state vector of the sub-harmonic mode in a 
squeezed number-basis 

|T ? ,n)=S(t ? )|n) J (25) 

with a dynamically adapted parameter rj = rj(t). 

The use of a dynamically optimized squeezed and dis- 
placed number basis [fj"2| allowed a numerical integration 
of the two-mode Schrodinger equation for input photon 
numbers up to several thousands. In Fig. 3 we have 
shown the scaled real part of the pump mode amplitude 
((02) = X2 +1P2) and its fluctuations as a function of the 
scaled time r = K (0,2(0)) t. Also shown is the mean-field 
result. One recognizes good agreement of the predictions 
for the optimum conversion time from both approaches. 




FIG. 3. Scaled in-phase average quadrature component of 
pump mode (S2) ~ (^2) / '(0:2(0)) from mean-field approxima- 
tion (line) and numerical integration of two-mode Schrodinger 
equation (long dashes). Also shown are the fluctuations 
of X2 obtained from numerical integration (short dashes), 
r = K(a 2 (0))t, (02(0)) = V200 

On the other hand, the numerical solution shows, that 
at the point of vanishing coherent amplitude of the pump 
mode, its fluctuations become large. This implies that 
the coherent-state approximation used in the mean-field 
approach is not valid near the point of maximum con- 
version. Furthermore, although the coherent amplitude 
vanishes, the mean photon number of the pump remains 
finite and thus the conversion efficiency is less than unity. 
Fig. 4 illustrates this. 
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fluctuations in the pump mode via the nonlinear interac- 
tion. As a result a finite amount of energy remains in this 
mode even at the point of vanishing coherent amplitude. 
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FIG. 4. Time evolution of energy of pump (dashed) 
and sub-harmonic mode (line) as function of scaled time 
r = K{a 2 (0))t for (o 2 (0)) = V200. 

Shown are the energies of both fields as a function 
of time. One recognizes a maximum conversion of only 
about 65% for 200 input photons of the pump mode. 
Our calculations indicate that this value does not increase 
with increasing input photon number and is thus not a 
finite-size effect. Near the point of maximum conversion 
the pump-mode amplitude becomes small and the back- 
action of the quantum fluctuations of the sub-harmonic 
mode (more precisely that of the anti-squeezed quadra- 
ture component) onto the pump mode gain importance. 
They lead to an increase of the in-phase quadrature fluc- 
tuations of the pump field and thus a finite amount of 
energy remains in this mode even though the coherent 
amplitude vanishes. 
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V. SUMMARY 



We have analysed the long-time quantum dynamics 
of degenerate parametric down conversion, for which 
standard approaches like the linearization of the Heisen- 
berg equations of motion fail. In a mean-field approach, 
which assumes a coherent pump mode but takes the sub- 
harmonic Auctions fully into account, an oscillatory en- 
ergy exchanges between the modes is found. The mean- 
field approach allows to determine the optimum inter- 
action times or crystal lengths for maximum squeezing 
and maximum down conversion. Since this approach ne- 
glects the quantum fluctuations of the pump, it becomes 
invalid near the point of maximum conversion and can- 
not be used to estimate the conversion efficiency. To 
calculate the latter we numerically integrated the two- 
mode Schrodinger equation. The numerical integration 
was possible for photon numbers up to several thousands 
due to the use of a dynamically optimized, displaced and 
squeezed number basis |l2j . We found that the maximum 
conversion efficiency is only about 65% for a coherent in- 
put of the pump mode. This limitation is a pure quantum 
effect. The large fluctuations in the anti-squeezed com- 
ponent of the sub-harmonic field introduce corresponding 
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